Kondo effect in quantum dots coupled to ferromagnetic leads 
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We study the Kondo effect in a quantum dot which is coupled to ferromagnetic leads and analyse 
its properties as a function of the spin polarization of the leads. Based on a scaling approach we 
predict that for parallel alignment of the magnetizations in the leads the strong-coupling limit of the 
Kondo effect is reached at a finite value of the magnetic field. Using an equation-of-motion technique 
we study nonlinear transport through the dot. For parallel alignment the zero-bias anomaly may be 
split even in the absence of an external magnetic field. For antiparallel spin alignment and symmetric 
coupling, the peak is split only in the presence of a magnetic field, but shows a characteristic 
asymmetry in amplitude and position. 

PACS numbers: PACS numbers: 75.20.Hr, 72.15.Qm, 72.25.-b, 73.23.Hk 



The Kondo effect Q in electron transport through a 
quantum dot (QD) with an odd number of electrons is 
experimentally well established 0, 0] • Screening of the 
dot spin due to the exchange coupling with lead elec- 
trons yields, at low temperatures, a Kondo resonance. 
The main goal of the present work is to investigate how 
ferromagnetic leads influence the Kondo effect. In the ex- 
treme case of half-metallic leads, minority-spin electrons 
are completely absent, i.e., the screening of the dot spin 
is not possible, and no Kondo-correlated state can form. 
What happens, however, for the generic case of partially 
spin polarized leads? How does the spin-asymmetry af- 
fect the Kondo effect? Is there still a strong coupling 
limit, and how are transport properties modified? 

Based on a poor man's scaling analysis we first show 
that the strong-coupling limit can still be reached in this 
case if an external magnetic field is applied. This is famil- 
iar from the Kondo effect in QDs with an even number 
of electrons 0, 0, El Q i which occurs at finite magnetic 
fields, although the physical mechanism is different in 
the present case. In the second part of the paper we an- 
alyze within an equation-of-motion (EOM) approach the 
nonlinear transport through the QD. We find that for 
parallel alignment of the lead magnetizations the zero- 
bias anomaly is split. This splitting can be removed by 
appropriately tuning the strength of an external mag- 
netic field B. In the antiparallel configuration of the lead 
magnetizations no splitting occurs at zero field. 

The Anderson Hamiltonian for a QD with a single level 
at energy ep coupled to ferromagnetic leads is 
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and in the QD, V r k is the tunneling amplitude, S z — 
(dtdi — d\d{)/2, and the last term is the Zeeman energy 
of the dot. (Stray fields from the leads are neglected.) We 
assume identical leads and symmetric coupling, Vj^k — 
Vru- The ferromagnetism of the leads is accounted for 
by different densities of states (DOS) v r ^(oj) and u r ±(uj) 
for up and down-spin electrons. 

In the following we study the two cases of parallel (P) 
and antiparallel (AP) alignment of the leads' magnetic 
moments. For the AP configuration and zero magnetic 
field and bias voltage, the model is equivalent (by canon- 
ical transformation |8() to a QD coupled to a single lead 
with DOS z/j,f + i/Rf = i>Li + Vri- In this case, the usual 
Kondo resonance forms, which is the same as for non- 
magnetic electrodes 

This changes for the P configuration. In this case, 
there is an overall asymmetry for up and down spins, say 
z/£f + j/Rf > vli + vri . To understand how this asymme- 
try affects the Kondo physics we apply the poor man's 
scaling technique 0, performed in two stages 0. In 
the first stage, charge fluctuations dominate and lead to 
a renormalization of the QD's levels. Since the renormal- 
ization for the spin-down level is stronger than for spin- 
up, a level splitting between the two spin orientations is 
generated. This is one of the key mechanism for all the ef- 
fects discussed below. To reach the strong-coupling limit 
it is, therefore, essential to apply an external magnetic 
field to compensate for the generated spin splitting. In 
the second stage, the resulting model is mapped onto a 
Kondo Hamiltonian and the degrees of freedom involving 
spin fluctuations are integrated out. For simplicity we as- 
sume for the scaling analysis flat DOSs v ra and neglect 
the fc-dependence of the tunnel amplitudes V r k = V. 



where c r ka and d a are the Fermi operators for electrons 
with wavevector k and spin a in the leads, r = L, R, 



First we reduce the cutoff D from Dq, which is the 
smaller value of the bandwidth and the onsite repulsion 
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U 01 • Charge fluctuations lead to the scaling equations 
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where we defined T a — 2ir\V\ 2 ^ r v rai and a is oppo- 
site to a. This yields the solution AT — Tf — Tj_ = 
-(l/7r)Prin(D /£ ) ) + Aeo, where P = (i/ rT - i/ rl )/(i/ r f + 
i/ r n measures the spin polarization in the leads, T — 
(T| + T|)/2 and Aeo = gfJ-BB is the Zeeman splitting. 
The empty-dot state |0) hybridizes with states where the 
dot is singly occupied |lcr) with either spin up or down, 
while the singly-occupied state \la) only hybridizes with 
the empty-dot state |0) (for U ^> |e| - asymmetric An- 
derson model). Because of the spin-dependent DOS in 
the leads the hybridization is spin-dependent, which is 
the physical origin of the generated AT. 

To describe Kondo physics (for T < 0) we terminate 
[lfij | the scaling of Eq. (J2J) at D ~ —e, and perform a 
Schrieffer- Wolff transformation. Using the renormalized 
parameters D and T we get the effective Kondo Hamilto- 
nian 
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plus the term — 2S z D(v^ J z ^ — i/j J z [) and a potential scat- 
tering term. The initial values for the coupling constants 
are J + = J_ = J z ^ = J z \_ = \V\ 2 /\e\ = J . To reach 
the strong-coupling limit we tune the external magnetic 
field B such that the total effective Zeeman splitting van- 
ishes, AT = (the field B will also slightly modify the 
DOS in the leads |(|). During the second stage of scal- 
ing spin fluctuations will renormalize the three coupling 
constants J + = J_ = J±, J z ^, and J z y differently. The 
scaling equations are 
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d{v±J±) 
d\n(D/D) 

d\n(D/D) 

with v± = yfv\ P| , v a = J2 r Vrcr E3- To solve these 
equations we observe that {v± J±) 2 — {i'iJzi)(t / iJzi) = 
and i/fj z f — v iJzi — JaP(y^ + v{) is constant as well. 
I.e., there is only one independent scaling equation. All 
coupling constants reach the stable strong-coupling fixed 
point J± — 3 Z \ — Jzi — oo at the Kondo energy scale, 
D ~ ksTx- For the P configuration the Kondo temper- 
ature in leading order, 
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depends on the polarization P in the leads. It is a max- 
imum for nonmagnetic leads, P = 0, and vanishes for 
P-yl. 



Finally, we point out an interesting consequence of 
the spin polarization in the leads. With nonmag- 
netic leads, the Kondo Hamiltonian couples the spin 
of the QD to the spin of the leads only, but not 
to its charge. To analyze the analogous situation 
in our case, we introduce the (pseudo) spin a = 

(V 2 )Efcfc'rrW clra<r™' c k>r>a>/(2Dy/vM), where the 
spin-dependent normalization factor is crucial to ensure 
the proper spin commutation relations, and the (pseudo) 
charge en = e J2kk'rr'a c kra c k'r'a/ (2Dv a ). The last term 
in Eq. (|3J can, then, be written as 2D(y^ J&\-\-v^ Jzl)&zSz 
plus D(ufJzf — v\J z AriS z . The first term is analogous 
to the Kondo model with nonmagnetic leads, while the 
second term couples spin to charge. The latter does not 
scale up and the associated additional renormalization 
of the Zeeman splitting, — (I/vr)Pr, is negligible as com- 
pared to AT in the limit Dq 3> |e|. 

The unitary limit for the P configuration can be 
achieved by tuning the magnetic field appropriately, as 
discussed above. In this case, the maximum conductance 
through the QD is G^ ax a — e 2 jh per spin, i.e., the same 
as for nonmagnetic leads. This yields that the ampli- 
tude of the Kondo resonance for up and down spin at the 
Fermi level are different, since G^ ax a ~ F ct (Pf)p<t(Pf), 
and therefore, p t (P f )M(Pf) = (1-P)/(1+P). For the 
AP configuration, the maximal conductance is reduced, 
GmSc.ff = (1 _ P 2 )e 2 /h, and vanishes for P — ► 1. 

In the remainder of this article we analyze the DOS 
of the QD and address nonequilibrium transport. For 
a qualitative discussion we should employ the sim- 
plest technique which accounts for both the formation 
of Kondo resonances and the influence of the spin- 
dependent renormalization of the dot level on spin fluc- 
tuations. The equation-of-motion (EOM) technique with 
the usual decoupling procedure |T3 . for higher order 
Green functions (GF) satisfies the first requirement but 
not the second. We, therefore, extend this scheme by 
calculating the level splitting AT self-consistently. For a 
more quantitative analysis one could include higher-order 
(than usually) Green's functions in the EOM approach 
or higher-order diagrams in the resonant-tunncling- 
approximationllj] , or use more advanced schemes such 
as real-time |l7j or numerical RG methods. These 
techniques are, however, much more complex |l5| and are 
not pursued here. 

Within the Keldysh formalism, the transport current 
I = J2 a la through a QD for T Ra (u) = A CT r L( r(w) is 

h j r LCT (w) + r R(7 (a;j 

where po-(w) = — (1/tt) ImG^ et (w). For strong interaction 
(U — > oo) the retarded Green's function is 



I - (ng) 



(8) 
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where E 0ct (w) = J^keL.R \Vk\ 2 / ~ £ ka) is the self-energy 
for a noninteracting QD, while 
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(9) 



appears for interacting QDs only. The average occupa- 
tion of the QD with spin a is obtained from (n a ) = 
— (i/2tt) J duo G< (lu) . Extending the standard derivation 
[12) , we replaced on the r.h.s. of Eq. © Ae — * A?, where 
e a is found sclf-consistently from the relation 



( a = e<x + Re[S 0CT (e C r) + Si^e^, Ae)] 



(10) 



which describes the renormalized dot-level energy, where 
the real part of the denominator of Eq. (JHJ) vanishes . 
We emphasize that without this self-consistency relation 
the Kondo resonances will, in general, appear at different 
positions, which disagree with the conclusions from the 
scaling analysis (lif . The procedure simulates higher- 
order contributions and the influence of the renormal- 
ization of the dot level on spin fluctuations. Follow- 
ing Ref. |l2l ] we introduce, in heuristic way, a lifetime 
Tct(/*l, MRj CT'^l) which describes decoherence due to a 
finite bias voltage V or level splitting Ae. It is obtained 
in second-order perturbation theory and depends on the 
electrochemical potentials in the leads, //l(r)j an d the 
level positions. Again, we replace the bare levels by the 
renormalized ones. In the numerical results presented 
below we use Lorentzian bands of width D = lOOR 

For nonmagnetic leads, P = and zero magnetic field, 
-6 = 0, the proposed approximation is identical to the 
standard EOM scheme 12]. For finite magnetic field, 
B 0, the self-consistency condition yields a splitting 
of the Kondo resonances which is slightly smaller than 
2<7/zb-B, in agreement with both experimental 0] and the- 
oretical findings [U HJ. For B = and P > in the 
parallel configuration, we obtain a value of the splitting 
A? comparable to the result from scaling, Eq. . 

In Fig. □ we plot the DOS of the QD for spins in AP 
and P configurations with spin polarization P = 0.2 in 
the leads. In the AP configuration there is one Kondo res- 
onance [Fig.QJa)] and the DOS is the same as for the case 
of nonmagnetic leads. For the P configuration, however, 
the Kondo resonance splits [Fig. [Jc)] , which can be com- 
pensated by an external magnetic field B [Fig. dd)]. In 
the latter case, the amplitude of the Kondo resonance for 
spin down significantly exceeds that for spin up (as dis- 
cussed above). A finite bias voltage, /xr, — = e,V > 0, 
again leads to a splitting for both the AP and the P con- 
figurations, [Fig.^b) and (e)]. In the AP configuration, 
the amplitude of the upper and the lower Kondo peak 
appear asymmetric [Fig.^b)]. 

In Fig. |3 we show the differential conductance as a 
function of the transport voltage. For nonmagnetic leads, 
there is a pronounced zero-bias maximum [Fig. Ufa)], 
which splits in the presence of a magnetic field [Fig.^b)]. 



For magnetic leads and parallel alignment, we find a 
splitting of the peak in the absence of a magnetic field 
[Fig. Etc)], which can be tuned away by an appropriate 
magnetic field [Fig. Efd)]. In the AP configuration, the 
opposite happens, no splitting at B = [Fig. Et e )] but 
finite splitting at B > [Fig. E£f)] with an additional 
asymmetry in the peak amplitudes as function of the bias 
voltage. This asymmetry is related to the asymmetry in 
the amplitude of DOS [Fig. Gib)] . All these findings are in 
good agreement with our scaling analysis. In Fig.[5{g) we 
show the tunnel magnetoresistance (TMR) which can be 
much larger than for conventional TMR systems. Finally, 
we find that the positions of the peaks in the AP con- 
figuration in the presence of a magnetic field are slightly 
shifted as a function of the polarization P [Fig. El n 
This can be explained in the similar way as in Ref. 
by an additional level splitting SAe — (l/47r) PT/eeV at 
finite bias voltages due to spin accumulation in the QD. 

We finally comment on the observability of our pro- 
posal. Ferromagnetic systems usually have strong spin- 
orbit coupling which could lead to spin-flip relaxation 
and, thus, could destroy the Kondo effect. On the other 
hand, recent observation of Aharonov-Bohm oscillations 
in a ferromagnetic ring |22) proved coherent transport, 
i.e., the formation of the Kondo cloud in ferromagnetic 
leads might be possible as well. How can one attach 
ferromagnetic leads to a QD? A conceivable realization 
might be to put carbon nanotubes in contact to ferro- 
magnetic leads jl^ ■ The Kondo effect has been observed 
already for nonmagnetic electrodes 0,13] • Alternatively 
one might use magnetic tunnel junctions with magnetic 
impurities in the barrier, or spin-polarized STM |25tl26| . 

In conclusion, we presented a qualitative study of the 
Kondo effect in QDs coupled to ferromagnetic leads. In 
particular we found a splitting of the Kondo resonance 
for parallel alignment of the leads magnetizations, even in 
the absence of a magnetic field. Our results are based on 
a poor man's scaling approach and an EOM technique. 
Further investigations on a more quantitative level using 
more advanced techniques would be desirable. 
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